Current Carrying States in a Random Magnetic Field 
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We report results of a numerical study of noninteracting electrons moving in two dimensions, in 
the presence of a random potential and a random magnetic field for a sequence of finite sizes, using 
topological properties of the wave functions to identify extended states. Our results are consistent 
with the existence of a second order localization-delocalization transition driven by the random 
potential. The critical randomness strength and localization length exponent are estimated via a 
finite size scaling analysis. 
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The problem of Anderson localization in two dimen- 
sions (2D) continues to be the subject of extensive re- 
search recently. Scaling argumentso predict that all one- 
particle states are localized in 2D, in the presence of scat- 
tering potential which respects time reversal symmetry. 
In the presence of a uniform magnetic field, however, be- 
cause of the fact that the magnetic field breaks time re- 
versal symmetry and suppresses the enhanced back scat- 
tering, extended states are found at individual critical 
energies!] These extended states are crucial to the exis- 
tence of the integer quantum Hall effect .a 

Recently there has been considerable interest in the 
localization problem in 2D in the presence of a random 
magnetic field with zero mean, motivated by the wor, 
of Halperin, Lee and Read,cl and Kalmeyer and Zhang J 
on the problem of the half filled Landau level. These 
authors argue that if one treats the electrons as compos- 
ite fermionsO carrying two flux quanta, the flux carried 
the the composite fermions cancels that of the external 
magnetic field on average when the Landau level filling 
factor is 1/2, and the systems may be viewed as a Fermi 
liquid like metal. Disorder induces local fluctuations in 
the density of the composite fermion, and hence at the 
mean field level the composite fermions see a fluctuating 
local magnetic field with zero mean, in the presence of 
disorder. 

Following this development, there have been many ef- 
forts devoted to the study of the localization problem 
of noninteracting particleS|-.rac)ving in a random mag- 
netic field with zero meanJ3~li3 However, the key issue, 
namely whether extended states exist in the thermodi^ 
namic limit, remains unresolved. Zhang and ArovastTJ 
argued that in the field theory description of the prob- 
lem, the local magnetic flux gives rise to a local topo- 
logical density which has long range interaction and in- 
duces a localization-delocalization transition (see, how- 
ever, Ref. [ [n]). Their conclu^n_is-.supported by nu- 
merical work of several groupajHOii3 which found ex- 
tended states near the band center on lattice models with 
weak random potential and finite gandom magnetic field, 
although Sugiyama and Nagaosafl concluded otherwise. 
On the other hand, Lee and ChalkeJa simulated a net- 
work model which they argue to be appropriate to the 
problem of random magnetic field, and find states arc lo- 



calised at all energies. Subsequently Kim, Furusaki and 
Leai-3 showed that the network model of Lee and Chalker 
may be mapped onto an SU(2N) spin chain in the limit 
N — > 0, which has short range correlations and an exci- 
tation gap; they thus confirm the conclusion of Lee and 
Chalker within the framework of the network model. The 
network model was originally introducedEJ in the study of 
the integer quantum Hall effect, and has enjoyed tremen- 
dous success in that context. Such success is ensured by 
the fact that in the presence of a strong uniform magnetic 
field, the magnetic length is very small, thus one body 
states live along one dimensional equipotential lines. In 
the case of a random magnetic field, however, the relevant 
percolating equal field lines are along B = 0, where the 
magnetic length is infinite] Therefore it is not completely 
obvious that the network model of Lee and Chalker is 
an appropriate description of the random magnetic field 
problem. 

In this paper we report a systematic finite size nu- 
merical study of the electron localization problem in the 
presence of a random magnetic field with zero mean, as 
well as a random potential, on the square lattice. We 
study topological properties of one-electron wave func- 
tions by calculating their, topological quantum numbers 
called the Chern number .1313 This approach has proved 
very successful in thp_stujdy of integer quantum Hall ef- 
fect and transitions E3~t3 In particular, finite size cor- 
rections appear to be much less in this approach, com- 
pared with others, implying a rapid convergence to the 
thermodynamic limit. Recently it has been applied by 
Sheng and Wengo to the present problem. There are, 
however, several important differences between our work 
and theirs.,-] ELcatly, their work as well as other numeri- 
cal studies Ja~l3'E2l concentrate on systems with a relatively 
weak (or zero) random potential, and attempt to identify 
a mobility edge. When the strength of random potential 
is weak, the length scales involved may be quite large, and 
it is often hard to distinguish whether the localization 
length near the band center is truly infinite or very large 
but finite. In contrast, we study systems with a series 
of different random potential strengths, with the same 
strength of random flux. The idea is that if extended 
states do exist for weak randomness, there must exist a 
critical random potential strength (which is typically not 
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very weak) at which all extended states disappear, and 
scaling behavior should be observable near this critical 
point. Secondly, unlike Sheng and Weng, who study be- 
havior of density of states (p(E)) and density of current 
carrying states (p c (E)), which are functions of energy E, 
we focus on quantities like the total number of extended 
states N c in a finite size system. iV c is the zeroth mo- 
ment of the distribution function p c {E). In other random 
systems such as spin glasses low order moments of distri- 
bution functions are known to converge to an acceptable 
level of accuracy relatively fast with sample averaging; 
the higher moments and the distributions themselves are 
still quite noisy at that stage. Finally, our massively par- 
allel computer allows us to study samples of considerably 
larger sizes than those in Ref . [ |l^] . 

We consider the one body tight binding Hamiltonian 
on the square lattice: 

H = — 5>^ct Cj + e-^cjci) + J2 e * c K (1) 
(y) « 

where Cj is the fermion operator on lattice site i. The 
first term represents hopping or kinetic energy of the elec- 
trons, and the summation is over nearest neighbors. The 
flux through each square is equal to the summation of 
Oij along its four edges. We take the flux to be ran- 
dom and uniformly distributed from —hir and fin, where 
< h < 1. The second term represents a random onsite 
potential. We take ej to be uniformly distributed between 
—w and w. h and w parameterize the strength of random 
magnetic field and random potential respectively; when 
h = 1 the flux through each square is completely random. 
We have set the hopping matrix element to be 1 and use 
it as the unit of energy. We study samples of square ge- 
ometry with linear size L (with number of sites N s = L 2 ), 
and impose periodic boundary conditions in both direc- 
tions: $(k + Lx) = e^^O), and ^(k + Ly) = e**(fc). 
The Hall conductance of an individual eigenstate \m) can 
be obtained easily using the Kubo formula: 

m _ ie 2 h (m\v y \n}(n\v x \m) - (m\v x \n) {n\v y \m) 
a *y- Ns (E n -E m y 

n^tra 

(2) 

where E n is the energy of the nth state, and 

v T = (i/h) 5>j +T c,e^+- - ctc i+T e-^+-') (3) 

3 

is the velocity operator along direction r (= x or y). 
a xy is identically zero in the absence of magnetic flux, 
as guaranteed by time reversal symmetry. In the pres- 
ence of random magnetic flux with zero average, the time 
reversal symmetry is broken in a specific configuration 
of randomness, and individual states may have nonzero 
o~ xy , although the disorder averaged Hall conductance is 
always zero for any given Fermi energy, since the aver- 
aging process restores the time reversal symmetry, tr™ 



depends on the boundary conditions imposed. As shown 
by Niu et at, the boundary condition averaged Hall con- 
ductance takes the formic 

«) = 4^2 / #1#2<(<^1, <h) = C(m)e 2 /h, (4) 

where C(m) is an integer called the Chern number of 
the state |m). States with nonzero Chern numbers carry 
Hall current and are necessarily extended statesB9 Thus 
by numerically diagonalizing the Hamiltonian on a grid 
of tf>i and (f>2 , and calculating the Chern numbers of states 
of finite size systems by converting the integral in (2) to 
a sum over grid points, we are able to identify extended 
states in a finite size system unambiguously. 

In this work we have studied finite size samples with 
linear size ranging from L — 4 to L — 19, for differ- 
ent random potential and random field strength. For 
each particular randomness strength and system size, 
we average over 30 to 2000 different configurations of 
randomness. Depending on system size and randomness 
strength, the number of grids necessary to determine the 
Chern numbers unambiguously using Eq. (^J) ranges from 
35 x 35 to 70 x 70. 

We plot in Fig. 1 the dependence of disorder averaged 
number of current carrying states N c versus the system 
size N s , for h = 0.5 and a series of different w, on a double 
logarithmic plot. We find for a given system size, N c de- 
creases as w increases. This is exactly what is expected, 
since random potential tends to localize states. For very 
strong randomness (w > 4.5) and big enough system size, 
we see iV c decreases with increasing N s , indicating there 
are no extended states in the thermodynamic limit, and 
the localization length is finite throughout the band. For 
somewhat weaker randomness, e.g. w = 4.0, although 
N c increases with N s within the size range we are able 
to reach, the curves obviously have downward curvature, 
and beyond certain N s , N c presumably decreases with 
increasing size and eventually goes to zero. Current car- 
rying states appear in finite size systems in this case be- 
cause the localization length £ for states near the band 
center, although finite, may be larger than the system size 
L. El For even weaker randomness (w < 3.0), however, we 
do not see any evidence of such downward curvature, and 
it appears that N c — > oo as N s — > oo. Such behavior is 
expected if extended states exist in the thermodynamic 
limit, either at individual critical energies, Elc3 or in a 
finite width band near the band center. In both cases, 
one expects the asymptotic behavior on the double loga- 
rithmic plot tp be linear. In the former case, the slope is 
expected to be less than 1, while for a band of extended 
states we have -/V c cx N s , hence the slope equals unity. 
For weak enough random potential, w < 2.0, we do find 
that the slope y is extremely close to 1, consistent with 
a band of extended states-in the thermodynamic limit. 

It has been suggestedllil that for weak randomness 
states are just barely localized (as in the case without 
magnetic fielda) with exponentially large £ near the band 
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center, much larger than system sizes accessible to nu- 
merical studies. While this may give the impression of 
an extended band, one would expect a slow, logarithmic 
decrease in N c /N s with increasing system size. However, 
we find no such evidence for small w; N c /N s appears to 
be a constant, or increases slightly. If, on the other hand, 
there is a band of extended states for small w, the width 
of such a band should decrease as w increases and shrink 
to zero at a critical random potential w c , above which £ 
is finite throughout the band (see Fig. 2). For w < w c 
we should have N c oc N s as N s — > oo, while for w > w c , 
N c — ► in this limit. Right at the critical point w = w c , 
there is a single critical energy at which the localization 
length .diverges , just like the case of a uniform magnetic 
field and we should have N c ~ N^", where y c is an 
exponent, related to the localization length exponent v 
throughErB 
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For w close to w c , there is a characteristic length scale 
that scales as £ m ~ \w — w c \~ u . For w > w c , £ m is 
the maximum localization length of the system, while for 
w < w c , £ m is the length scale above which the behavior 
N c oc N s is seen. Near the critical point we should see 
scaling behavior: 



N^F ± (L/U) ~ N^F ± (N s \w - w c \ 2 n, (6) 



where F^(x) are two different scaling functions for w < 
w c and w > w c respectively, which should satisfy the 
asymptotic behavior F + (x) ~ x 1 ~ Vr - and F~(x) — > as 
x — > oo. Therefore, if there is indeed a critical random- 
ness w c , we ought to see scaling behavior (^J), by tuning a 
single parameter w c . The other parameter v in the scal- 
ing relation (||) may be determined using Eq. (||), since 
the slope y is know from the data presented in Fig. 1. 

Such scaling behavior is indeed seen, and the data col- 
lapses best for w c = 3.0 and v = 1.8, as shown in Fig. 
3. At w = 3.0 the linear fit of log Y c versus log Y s yields 
y = 0.87±0.01, and from Eq. §) we obtain v = 1.9±0.2, 
consistent with scaling results. We emphasize that data 
collapsing of two different scaling curves is achieved by 
tuning one parameter, due to the constraint Eq. (|^). 
The scaling behavior we see here supports the existence 
of a band of extended states for weak random potential, 
and would not be expected if there were no localization- 
delocalization transition at finite w c . In order to make 
sure that our data are indeed in the scaling regime, we 
plot in Fig. 4 both the total density of states p(E) and 
density of current carrying states p c {E) (per site), for 
w = 3.0 and h = 0.5, for two different system sizes L = 4 
and L = 16. We find the width of p c (E) is considerably 
smaller than that of p(E), especially for the larger size. 
Further the width of p c (E) is size dependent while that of 
p(E) is essentially the same for both sizes. This is similar 
to what is seen in the-s.tudy of the localization problem in 
lowest Landau level,Eil and indicates that we are indeed 



in the scaling regime. By contrast the width of p c (E) 
and p{E) is almost the same for Gaussian white noise 
potential in the first Landau level, W^Gfe the localization 
lengths are known to be very largeBcJ 

From the scaling we are also able to estimate the crit- 
ical randomness to be w c = 3.0 ± 0.2, and v = 1.8 ± 0.2, 
for h — 0.5. The localization length exponent v we ob- 
tain bore is very close to that estimated by Sheng and 
Wenglj using a very different scaling scheme, although 
they do not have a quantitative estimate of the critical 
random potential strength w c . 

In summary, we have found good evidence supporting 
the existence of a localization-delocalization transition 
driven by random potential, in the 2D square lattice in 
the presence of a random magnetic field with zero aver- 
age, implying the existence of extended states at weak 
randomness. The critical random potential strength and 
localization length exponent are estimated using finite 
size scaling analysis. 
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N s , for h = 0.5 and various w on a double logarithmic scale. 
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FIG. 3. The scaling functions F ± (N a \w - w c \ 2 "). 

FIG. 4. Disorder averaged density of states p(E) and den- 
sity of current carrying states p c (E) at w = 3.0 and h = 0.5, 
for L = 4 and L = 16. 
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